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Abstract We present a system for molecular spec- 
troscopy using a broadband mid-infrared laser with near 
infrared detection. Difference frequency generation of a 
Yhffbcr femtosecond laser produced a mid-infrared (MIR) 
source tunable from 2100 - 3700 cm" 1 (2.7 - 4.7 /im) 
with average power up to 40 mW. The MIR spectrum 
was upconverted to near-infrared wavelengths for broad- 
band detection using a two-dimensional dispersion imag- 
ing technique. Absorption measurements were performed 
over bandwidths of 240 cm" 1 (7.2 THz) with 0.048 cm" 1 
(1.4 GHz) resolution, and absolute frequency scale un- 
certainty was better than 0.005 cm -1 (150 MHz). The 
minimum detectable absorption coefficient per spectral 
element was determined to be 4.4 x 10~ 7 cm -1 from 
measurements in low pressure CH4, leading to a de- 
tection limit of 2 parts-per-billion. The spectral range, 
resolution, and frequency accuracy of this system show 
promise for determination of trace concentrations in gas 
mixtures containing both narrow and broad overlapping 
spectral features, and we demonstrate this in measure- 
ments of air and solvent samples. 



1 Introduction 

Molecular spectroscopy has long been a useful tool in 
chemistry and physics. As applications for molecular spec- 
troscopy continue to grow into environmental sensing [1 , 
security [2J, and medical analysis [3J, the requirements 
for broad spectral range combined with good sensitivity 
become more demanding. Traditional molecular spec- 
troscopy has been accomplished with thermal sources 
in moving mirror Fourier transform spectrometers or 
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with tunable single frequency lasers, but recent devel- 
opments in femtosecond laser frequency combs combine 
the high brightness and frequency precision of a tunable 
laser with the broad spectrum previously open only to 
thermal sources. 

Effective molecular spectroscopy requires both the 
source at a useful wavelength range as well as a de- 
tection method to utilize the unique advantages of the 
source. Broadband laser sources spanning 1.0 — 1.8 /im 
wavelengths are readily available but only able to ob- 
serve weaker overtones of the strong hydrogen-stretch 
absorption features of many target molecules. Sources in 
the mid-infrared wavelengths from 2.0 — 3.6 /urn would 
access much stronger fundamental absorption to allow 
more sensitive detection of these molecules. Broadband 
femtosecond MIR lasers have been demonstrated using 
difference frequency generation between a comb and CW 
laser [J], optical parametric oscillation [SJH5], as well as 
direct comb generation [TJ. 

Detection techniques that take advantage of the wide 
spectral bandwidth of a femtosecond source include us- 
ing the laser as a high brightness source for a moving 
mirror Fourier transform spectrometer [7J[S] or dual fre- 
quency comb Fourier transform spectroscopy [5 UT0"Illl| . 
Dispersive methods for imaging a frequency comb spec- 
trum have been shown with Vernier coupling of comb 
and cavity [12] or a high resolution etalon [3l H3lll4| . 
Matching the mode spacing of a cavity to the repetition 
rate of a mode-locked laser allows broadband cavity en- 
hanced absorption [15] or cavity ringdown spectroscopy 

H5I- 

We demonstrate a femtosecond laser system that gen- 
erates broad spectra covering an important mid-infrared 
absorption range of 2.7— 4.7 /im (2100— 3700 cm -1 ) while 
retaining advantages of detection in the near infrared 
through an upconversion process. A single measurement 
with up to 5000 elements at 0.048 cm _1 (1.4 GHz) reso- 
lution enables simultaneous measurement of sharp spec- 
tral lines and broad absorption features from complex 
molecules. 



2 



T. A. Johnson, S. A. Diddams 



1040nm MIR 800nm 




Fig. 1 A Yb:fiber femtosecond laser was separated, broad- 
ened, and recombined on a nonlinear crystal to produce broad 
bandwidth mid-infrared spectra via difference frequency gen- 
eration. The mid-infrared beam was passed through a gas 
cell and upconverted on a second nonlinear crystal, then 
transferred to the two-dimensional spectrometer for anal- 
ysis. PBS=polarizing beam splitter, DBS=dichroic beam- 
splitter, MF=microstructured fiber, SMF=single mode fiber, 
PPLN=periodically poled lithium niobate crystal. 

2 Experimental Details 

2.1 Broad bandwidth mid-infrared source 

A Yb:fiber femtosecond laser was the source for a dif- 
ference frequency generation (DFG) approach to a tun- 
able, broad bandwidth MIR source. The experimental 
setup is shown in Figure [1] Full details and character- 
ization of this source are provided in a separate publi- 
cation [17] , while here we provide a brief summary of 
the properties most relevant for our spectroscopic mea- 
surements. The MIR DFG portion of the system has 
similarities to Enfiber comb approaches [T51I19II2T)] . The 
output of an amplified Yb:fibcr femtosecond laser pro- 
vided 125 fs pulses at 100 MHz and 2.5 W average power. 
This source was split into pump and signal paths, seed- 
ing a fan-out, periodically-poled MgO : LiNbC>3 (PPLN) 
nonlinear crystal in a single pass to generate an idler 
wavelength (A" 1 - A^ 1 = \~r l ). The pump path deliv- 
ered 2.0-2.4 W average power to the PPLN after passing 
through a variable path delay stage. The signal path 
coupled 100-500 mW average power into aim length 
microstructured fiber to produce a Raman-shifted soli- 
ton that was tunable from 1.3 — 1.7 /im by varying cou- 
pled power (see Figure [2ji). The pump and signal were 
spatially and temporally overlapped in the PPLN to pro- 
duce the MIR light. Figure [^p shows the broad tuning 
possible by adjusting the in-coupled power to the mi- 
crostructured fiber and the PPLN phase matching po- 
sition. After the PPLN, the MIR was separated from 
the pump and signal beams using a broadband anti- 
reflection coated germanium filter. MIR generation ef- 
ficiency varied with pulse parameters and crystal wave- 
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Fig. 2 (a) Optical spectrum analyzer measurements of the 
filtered Raman solitons for average input powers 100-500 mW 
(average output power 10-20 mW). Each spectrum normal- 
ized to peak value, (b) Monochromator and InSb photodetec- 
tor measurements of the MIR source. Each spectrum normal- 
ized to peak value, with some spectra displaying absorption 
signals due to atmosphere (such as the carbon dioxide ab- 
sorption near 4.2 fim). Color of each signal spectrum in (a) 
corresponds to color of generated MIR spectrum in (b). 

length properties, and at the time of this work maxi- 
mum power up to 40 mW occured near A; ~ 3.3 /im. 
Observations reported in reference |T7] indicate coher- 
ence degradation of the Raman-shifted soliton and MIR 
light, however, the source had relatively low intensity 
noise (< 1%, 10 Hz - 10 MHz), a broad spectrum, and a 
well-collimated beam that were necessary for this work. 



2.2 Absorption path and gas samples 

The MIR beam was mode-matched into a multi-pass 
astigmatic Herriott cell (Figurc[TJ. The total path length 
was 210 m, and mirror losses on the 238 reflections re- 
sulted in an optical transmission of 2% for the cell. A 
gas handling system was attached to the cell to prepare 
different samples, with total cell and system volume of 
6 L. Spectral artifacts (due to mirror or coating spec- 
tral dependence in the beam path) were minimized by 
comparing the measured gas absorption features with a 
normalization measurement of the evacuated cell. 
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2.3 Virtually Imaged Phase Array (VIP A) and 
upconversion detection 

Rapid detection of broadband absorption signals was ac- 
complished with a two-dimensional dispersion method 
using a Virtually Imaged Phase Array (VIPA) crossed 
with a conventional diffraction grating [2"Tll2"2"] . The VIPA 
acted as a highly dispersive transmission grating with 

1.4 GHz (0.048 cm -1 ) resolution and free spectral range 
of order 100 GHz (3.3 cm -1 ). The overlapping orders 
from the VIPA were spatially separated by the crossed 
grating and recorded on a CCD camera (FiguresQ]and[Si) 
to analyze a large source bandwidth simultaneously. Spec- 
troscopic measurements with a VIPA have been demon- 
strated in the visible [13] and near-infrared [14] wave- 
lengths, but extension of the VIPA approach directly to 
the MIR has additional complications including optical 
material manufacturing limitations and the availability 
and sensitivity of cameras. 

To enable use of the existing near infrared VIPA 
and cooled silicon CCD camera for MIR spectroscopy, 
our system utilized an additional upconversion step to 
transfer the longer wavelength spectral absorption fea- 
tures to a detectable wavelength range [231121] . The MIR 
beam exited the gas cell with absorption features im- 
printed on the source spectrum and was combined on 
a dichroic beamsplitter with a CW 500 mW single fre- 
quency Nd:YAG laser at 1064nm (Figure]!]). Both beams 
were focused into a second PPLN crystal phase-matched 
for sum frequency generation, producing upconverted 
pulses near 800nm. The upconverted signal was trans- 
ferred in single mode fiber to the VIPA optics. The 
length of the upconversion PPLN was chosen to maxi- 
mize the available MIR and VIPA spectral ranges rather 
than total power, resulting in an average upconversion 
quantum efficiency of about 2 x 10~ 5 . Combining the gas 
cell mirror losses and upconversion efficiency with the 
additional near-IR optical losses in the transfer fiber and 
VIPA spectrometer, the silicon CCD camera recorded 
about one photoelectron per 10 9 MIR photons generated 
at the first PPLN crystal. Despite this significant loss in 
signal, it was still possible to recover useful spectral in- 
formation due to the advantage of strong MIR molecular 
absorption features rather than direct detection of weak 
overtones with a near-infrared source. 

2.4 Frequency calibration of VIPA images 

The frequency information of the absorption spectrum 
was spread across the CCD according to the dispersion 
law described in reference [25] , which accounts for VIPA 
thickness, tilt angle, and beam characteristics. In prac- 
tice it is difficult to assign frequencies directly from this 
dispersion law, and it is necessary to employ some form 
of direct image calibration by referencing known absorp- 
tion spectral features [33] or directly imaging the fre- 
quency comb source [12]. Calibration of the frequency 
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Fig. 3 Two-dimensional VIPA spectrometer images cap- 
tured by CCD. (a) Dispersed source spectrum with no gas 
sample in the cell, (b) Methane absorption features imprinted 
on source spectrum, (c) Subtracted image (a)-(b), to high- 
light the absorption spectrum pattern on the VIPA image. 
The VIPA's free spectral range is apparent in the vertical 
repetition of the absorption feature pattern at the top and 
bottom of image, (d) 10 GHz Ti:Saph comb calibration grid, 
optical filter etalon causes variations in comb mode intensity. 

information contained in the VIPA images was instead 
accomplished by launching a 10 GHz TkSaph frequency 
comb [2B] and CW 795 nm reference laser into the VIPA 
spectrometer along with the upconverted signal (Fig- 
ure[lj. When combined with measurements of the TkSaph 
repetition rate / rG p.TS and absolute frequencies of the 
CW reference (/79s) and upconversion pump (/1064) lasers, 
an absolutely calibrated frequency grid was correlated to 
the original upconverted VIPA image. The position of a 
particular mode of the TkSaph comb on the VIPA image 
(Figure ]3]i) was assigned a MIR frequency value 

/mode = /V95 + m x FSR + n x / rcp ,TS - fiou (1) 

where FSR was the free spectral range of the VIPA, m 
was the number of free spectral ranges (vertical stripes) 
separating the TkSaph comb mode and 795 nm laser po- 
sitions on the image, and n was the number of TkSaph 
comb modes counting up from the bottom of each free 
spectral range. Image positions between the TkSaph modes 
were assigned by linear interpolation between successive 
modes. The precise value of the VIPA FSR depended 
upon the spectrometer tilt angle, and was calibrated 
from the VIPA images based on comparison with / rep ,TS 
over the entire image range. The accuracy limit of / mo de 
was expected to be about 100 MHz due to the wave- 
length meter measurements of the CW lasers. In prin- 
ciple it would be possible to simplify frequency calibra- 
tion by directly resolving a MIR comb with either an 
increased MIR comb / rep or selective filtering of comb 
modes [13] in combination with improved VIPA optical 
quality. 



4 



T. A. Johnson, S. A. Diddams 



2. 5 Measurement procedure 

With transmitted light intensity I from incident inten- 
sity Jo through a sample length L, the Beer-Lambert 
Law predicts an absorption coefficient of 

To account for camera count offsets, optical losses, 
stray light, and variation of the source spectrum, a full 
measurement consisted of pumping out the cell to record 
an upconverted normalization spectrum (Figure^), fill- 
ing the cell and recording a VIPA image of the up- 
converted gas transmission spectrum (Figure [3b) , and 
blocking the MIR beam before the upconversion crystal 
to obtain a background signal level (not shown). Images 
of the reference CW lasers and Ti:Saph comb (Figure[3ji) 
were recorded to produce the frequency calibration grid, 
and each 2-D VIPA image was reconstructed into a set 
of frequency and camera count pairs. The absorption co- 
efficient a for each frequency step was calculated as 

1 , / Ctran — Cbg \ , , 

a = -J ln [r 7^1 (3) 

J- 1 \ Worm — <->bg / 

where Ctran, C norm , and Cb g are the counts in the trans- 
mission, normalization, and background images. 

3 Absorption spectra 

3.1 Low pressure CH4 

The evacuated Herriott cell was filled with natural iso- 
topic abundance methane gas (CH4) at 2.3 ± .1 mTorr 
total pressure. Figure S] shows the measured absorption 
spectrum and a predicted spectrum calculated from the 
HITRAN database [27] . The absorption spectrum was 
recorded at one MIR tuning position for the Raman 
soliton and DFG PPLN, using a single 5 s exposure 
for each of the transmission, normalization, and back- 
ground images. The 5 s exposure time was chosen to en- 
sure no CCD pixels exceeded half of the the maximum 
well depth (77,000 electrons). Total measurement cycle 
time of about 60 s was most limited by the pumping 
speed for gas exchange. 

At low pressure, linewidths were dominated by the 
optical resolution of the VIPA spectrometer (1.4 GHz or 
0.048 cm" 1 ), so the HITRAN-predicted spectrum was 
convolved with the optical profile of the instrument to di- 
rectly compare the spectra. The absolute frequency scale 
in this data contains no adjustable parameters. Com- 
parison of the measured and predicted frequency scales 
shows agreement within about 150 MHz (0.005 cm" 1 ), 
comparable to the expected uncertainty in frequency cal- 
ibration of / mo de- 
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Fig. 4 (a) Methane absorption signals (above axis) and pre- 
dicted absorption spectrum from HITRAN database (below 
axis) in 210 m path length and 2.3 mTorr pressure, (b) and 
(c) show narrow frequency regions of the the same data from 
(a). Linewidths were limited by optical resolution of the 
VIPA spectrometer, and the HITRAN spectrum was con- 
volved with the optical profile to compare directly to the 
measured spectrum. The entire spectrum was recorded at a 
single MIR tuning position in a single 60 s measurement cycle 
(normalization, transmission, and background images). 



Each VIPA order was recorded on the image (Fig- 
ure [3]) as a stripe with frequency increasing in the ver- 
tical direction. All pixels in a stripe at the same ver- 
tical position had the same frequency value and were 
integrated horizontally to provide the count number in 
a single frequency step. Each vertical pixel step on the 
image was separated by about 0.28 GHz (0.009 cm" 1 ), 
which was smaller than the optical resolution of 1.4 GHz 
(0.048 cm" 1 ). In the highest intensity portion of the nor- 
malization image, the integrated electron counts in each 
vertical step were C n orm ~ 65,000. The cooled CCD 
camera had low dark count and readout noise, so the 
dominant sources of noise were from shot noise in the 



Mid-infrared upconversion spectroscopy based on a Ybifiber femtosecond laser 



5 



electron count number (V C nor m and VCtran) and tech- 
nical noise due to drift of the laser intensity or optical 
alignment. Decreased signal to noise at the edges of the 
spectrum was due to decreased spectral intensity of the 
MIR source in these regions. 

The best signal to noise ratio for a single frequency 
step would be expected with transmission that is neither 
very high nor very low. At high absorption, the small re- 
maining counts Ctran would be indistinguishable from 
Cnorm noise, leading to a maximum detectable absorp- 
tion coefficient of a maXjS tcp ~ 2.6 x 10~ 4 cm -1 for our 
experimental conditions. At low absorption, the signal 
C norm — Ct ran would be comparable to the total noise 
from C norm and Ct ran , leading to a minimum detectable 
absorption coefficient of a m i niS tcp ~ 2.6 x 10~ 7 cm -1 . 
The range of a mi n, s tc P to a max , st cp represents an esti- 
mate based on a single frequency step at the highest 
intensity portion of the image. Changing the absorption 
path length L would alter the maximum and minimum 
observable absorption coefficients while maintaining the 
same dynamic range, and spreading the absorption infor- 
mation over additional pixels would allow an increased 
dynamic range by avoiding the pixel well depth limit. 

The limited VIPA optical resolution (1.4 GHz) re- 
sulted in a single spectral resolution element receiving 
contributions from about five 0.28 GHz frequency steps, 
so we would expect a statistical improvement of roughly 
l/y/5 on the noise level for each spectral resolution el- 
ement and a minimum detectable absorption coefficient 
of GVmin.ciom ~ 1-2 x 10 -7 cm -1 per element. In the case 
of low pressure methane, the natural linewidth was nar- 
rower than the spectral resolution, so the minimum de- 
tectable absorption coefficient for a single peak should 
correspond to approximately a m in,eiem- Additional tech- 
nical noise or analysis in lower intensity portions of the 
VIPA image would decrease the detectable absorption 
range. 

Another approach to determining the minimum de- 
tectable single peak absorption coefficient is based on the 
work of Hubaux and Vos [28) . In this technique, mea- 
surements of a single line in a set of different sample 
concentrations are used to determine the minimum de- 
tectable concentration level that is not consistent with 
zero. A modified version of this technique [14] uses the 
many available line strengths in a single recorded spec- 
trum to make this determination, and we employed a 
similar method to analyze our system. The predicted 
HITRAN spectrum for methane was prepared as a refer- 
ence, including the optical profile convolution. All locally 
prominent spectral features in the predicted spectrum 
with peak absorption coefficients > 10~ 7 cm -1 were se- 
lected as candidates. Some of these features were not sin- 
gle lines, but contributions from several closely spaced 
lines. 

Absorption coefficient baseline offset errors due to 
drift of the laser intensity or optical alignment would be 
particularly noticeable at high or low absorption coef- 



ficients, but measurement of isolated spectral lines al- 
lowed a method to partially account for this drift by 
using the regions between absorption features. Before 
fitting the absorption peak directly, a broader range of 
data around the peak (typically 4 cm -1 to 6 cm^ 1 ) was 
used to determine if a vertical offset was present in the 
measurement. The HITRAN model spectrum for this en- 
tire region was least squares fit to the data with a global 
amplitude adjustment to find the average baseline offset. 
While holding the offset value fixed and then using only 
the data within the peak's narrow width, the HITRAN 
model was fit to the peak a second time to find the mea- 
sured peak amplitude. The initially predicted pressure 
must be close to the actual pressure to limit significant 
changes in lineshape or position due to pressure broad- 
ening and shifting. The baseline offset of the absorption 
coefficient varied nearly linearly across the entire spec- 
trum by about 8 x 10~ 7 cm -1 , consistent with a spectral 
intensity drift of less than 2% between normalization and 
transmission images. 

The measurement fit peak amplitude was plotted ver- 
sus the peak amplitude of the original predicted HI- 
TRAN spectrum, and the process was repeated for the 
rest of the selected peaks. Figure [SJi shows the result- 
ing plot along with a linear fit. The slope of the linear 
fit finds the measured methane pressure compared to 
the predicted pressure, and the uncertainty in the slope 
provides an uncertainty in the measured pressure. The 
measured methane pressure was 2.233±.005 mTorr, con- 
sistent with the pressure gauge reading of 2.3± .1 mTorr. 
The 3<7 confidence intervals that correspond to this fit 
are also shown. The vertical intercept of the upper confi- 
dence interval can be interpreted as the minimum value 
of a peak's absorption coefficient that must be measured 
to have 99.87% confidence that the measured feature is 
not consistent with zero. The la vertical intercept (not 
shown) of a no ise ~ 4.4 x 10 -7 cm -1 corresponds to a 
peak height at the effective noise floor of the system (for 
a single spectral element). This noise level approaches 
within a factor of four of the shot noise limited detec- 
tion level ccmimeiem, and is equivalent to less than 1% 
fractional laser intensity noise. 

The 3<7 confidence intervals as presented here are for 
the measurement of a single peak within the spectrum. 
In the context of a concentration measurement, the sen- 
sitivity can be improved by using many different peaks 
to average down the noise on the individual peak de- 
tections. One implementation of the multi-line detection 
advantage is to include a statistical improvement in the 
Hubaux- Vos fit of 1/ y/N, where N is the number of spec- 
tral peaks within the measurement. This approach is evi- 
dent in [14] from the close proximity of the 3er confidence 
intervals to the data spread, although the resulting 3<r 
detection limit on the vertical intercept is more difficult 
to interpret as a practical figure of merit. Alternatively, 
comparison of the individual frequency step noise floor 
with the total integrated absorption in the measurement 
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Linear fit 3cr confidence interval 

(a) CH4 (b) H2O in 620 Torr air (c) CH4 in 620 Torr air 
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Fig. 5 Measured versus predicted absorption coefficients for gas samples. The linear fit slopes calibrate the measured pressure 
versus the predicted pressure, and the vertical intercept of the 3cr confidence interval shows the minimum measured absorption 
coefficient required to identify a single peak with 99.87% confidence that it is not consistent with zero. In all plots, the linear 
fits were constrained to pass through the origin. The total number of measured peaks is also shown for each graph. 



bandwidth [5] yields a theoretical concentration limit 
taking into account both the multi-line advantage as well 
as the varying absorption strengths of all lines. Using our 
measured system noise floor and bandwidth, this anal- 
ysis indicates a minimum concentration limit of about 
2 parts-per-billion (ppb) would be expected for methane 
in pure nitrogen at a pressure of one atmosphere. Addi- 
tional measurements would be necessary to verify this. 

3.2 Atmosphere 

The gas cell was vented to atmosphere and a sample of 
background air from the laboratory filled the cell to to- 
tal pressure 620 Torr. The recorded spectrum is shown 
in Figure \§\ with HITRAN predicted spectra for H2O 
(partial pressure 7.4 Torr) and CH4 (partial pressure 
1.1 mTorr), where the pressures for H2O and CH4 were 
initial estimates based on typical air composition and 
laboratory conditions. No significant absorption features 
were observed for other species. The entire span was 
again recorded at a single instrument position, with three 
5 s recorded images and about 60 s total acquisition time. 
Pressure broadening dominated the lincwidths over the 
optical profile. The measured spectrum began to show 
saturation as the absorption coefficient exceeded a sat 
1.8 x 10~ 4 cm -1 . As discussed in Section ETT1 this level 
of saturation could be attributed to fractional noise or 
drift of the comb intensity of order 2%. 

To limit the effect of the saturated peaks, any fre- 
quency step with a > 1.0 x 10~ 4 cm" 1 was removed 
from the analysis data before peak fitting. The offset 
and peak fitting protocol was similar to that described 
in Section 13.11 and allowed independent variation of the 
amplitude of the water vapor and methane signals in the 
fitted model. The fit was still performed in regions where 
saturated points were removed based only on the wings 
of the peak (provided that at least half of the points 
within the linewidth were still present). The measured 
versus predicted absorption coefficient comparisons with 



linear fits are shown in Figure [5}>c with partial pres- 
sures of 7.78 ± .06 Torr for H 2 and 1.27 ± .02 mTorr 
for CH4 that are within the expected range for the lab- 
oratory atmosphere conditions. The methane concentra- 
tion in a 620 Torr atmosphere can also be understood as 
2.05 ± .03 parts-per-million. 

The minimum detectable absorption coefficients for 
both methane and water vapor were found to be ten 
times higher than the low pressure methane measure- 
ment, despite identical apparatus and technique. This 
was also evident from the increased methane concentra- 
tion uncertainty of 30 ppb as compared to the predicted 
performance of about 2 ppb. This noise floor increase was 
due to a combination of effects related to a "real- world" 
gas sample that includes overlapped spectral features, 
saturated peak absorptions, and additional optical scat- 
tering from airborne particles. 



3.3 Solvent mixture 

A tradeoff exists between the small wavelength range 
with high signal to noise of a tunable CW laser versus 
the wide coverage of the broadband MIR/VIPA system 
with a decreased signal to noise for any single spectral 
line. The tunable CW laser would retain a detection ad- 
vantage when only a small number of isolated lines are 
present, but as the density of lines and species increase 
or overlap, it would become challenging to determine a 
concentration with a single CW laser. The broadband 
approach would benefit from additional lines within its 
span by using a collective fitting of all peaks to lower 
the detection limit for a gas species [8J. The difference 
between the CW laser and broadband laser becomes par- 
ticularly apparent when measuring a continuous absorp- 
tion spectra, and we show measurements of the solvents 
toluene (CyHg), gamma-butyrolactone (C4H6O2), and 
isophorone (C9H14O) to illustrate this possibility. These 
specific solvents were chosen because they have been 
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Fig. 6 Absorption spectrum of laboratory air at 620 Torr. The total measurement is shown above the axes, and HITRAN 
predicted spectra for methane and water vapor are shown below the axis. Saturation effects are evident for the strongly 
absorbing water vapor peaks. 



identified as some of the compounds commonly present 
in the headspace around explosive materials [29 . 

A gas sample mixed from all three solvents was pre- 
pared to study concentration measurements of broad- 
featured gases. The reference spectra for isophorone and 
gamma-butyrolactone were isolated from previous cali- 
bration with the MIR upconversion and VIPA spectrom- 
eter, and a high resolution quantitative absorption spec- 
trum was used for toluene |30] . The gas handling sys- 
tem was not designed to prepare this type of mixture 
precisely, so the partial pressures of each solvent were 
not well known. In addition, contamination of the liquid 
solvents by exposure to atmosphere during mixing intro- 
duced an unknown ratio of other gases to the mixture. 
Once expanded into the Herriott cell, the gas mixture 
had a total pressure of 232 ± 1 mTorr. 

The total absorption spectrum is shown in Figure [7J 
The narrow peaks relative to the broadly varying back- 
ground were compared to the HITRAN database and 
found to match H2O that was unintentionally introduced 
during mixing. No other narrow spectral features were 
found. The partial pressures of the three solvents and 
water vapor were calculated by applying a least squares 
fit to match the total measured spectrum. The indepen- 
dent absorption spectra of the constituent gases in the 
mixture are also shown in Figure [7] for the fitted par- 
tial pressures of 142 ± 6 mTorr (H 2 0), 18.6 ± .4 mTorr 
(toluene), 57±6 mTorr (isophorone), and 7.5±1.7 mTorr 
(gamma-butyrolactone). The largest sources of uncer- 
tainty in these pressures arose from the reference spec- 
tra used for calibration. The total pressure based on the 
spectral fits of the constituent gases was 225 ± 9 mTorr, 
which agrees with the total direct pressure measurement 
of 232 ± 1 mTorr. The subtraction residuals (Figure [7]) 
had a standard deviation of about 9 x 10~ 7 cm -1 for 
the individual frequency steps (not resolution elements), 
a value that is within a factor of four of the expected 
shot limited single frequency step value a m i n ,stcp ~ 2.6 x 




2900 2950 3000 

wavenumbers [cm 1 



Fig. 7 Measured total absorption of the solvent mixture 
(top). The other spectra show the contributions of the 
constituents found by fitting the total spectrum with ad- 
justable partial pressures. Absorption coefficient fit residuals 
are shown below the main graph and have a standard de- 
viation of about 9 x 10~ 7 cm -1 , within a factor of four of 
the shot noise limited minimum absorption coefficient for a 
single frequency step. 



10 7 cm 1 and corresponds to intensity noise levels of 
less than 2% on the normalization comb spectrum. 

The concentration precision in this demonstration 
was most limited by the non-ideal gas handling system 
for both the solvent mix and two of the reference spectra. 
A more fundamental limitation when measuring broad 
absorption features is lack of opportunity to compen- 
sate for drift through the baseline correction methods 
used in the low pressure CH4 and air measurements, re- 
quiring a more stable system on both the DFG process 
and upconversion detection process. Some of this drift 
could be eliminated by reducing the time interval be- 
tween the normalization and absorption images with a 
faster pumping speed. 
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4 Conclusion and outlook 

This system demonstrates the first high resolution spec- 
troscopic application of femtosecond lasers in the mid- 
infrared without the complication of Fourier transform 
approaches using moving optical paths or multiple combs. 
The 0.048 cm -1 (1.4 GHz) resolution spanned up to 
240 cm -1 (7.2 THz) in a single measurement to yield 
up to 5000 distinct spectral elements. The spectral el- 
ements exhibited a minimum detectable absorption co- 
efficient of a no i so « 4.4 x 10~ 7 cm -1 , a value within a 
factor of four of the expected shot noise limited absorp- 
tion coefficient. This level of absorption coefficient noise 
could be expected to approach a concentration detec- 
tion limit of 2 ppb for methane in pure nitrogen, and 
the present system has demonstrated methane concen- 
tration uncertainties at the level of 30 parts-per-billion 
in a more demanding direct atmospheric sample. 

There are clear paths to improving the detection sen- 
sitivity of this system. First, improved gas handling would 
increase measurement bandwidth by an order of magni- 
tude and likely reduce intensity drift noise between ab- 
sorption and normalization images. Second, it is possible 
to increase the MIR power up to 25 times [5] and the 
upconversion CW pump power more than 20 times to 
yield a 500-fold improvement in signal. Taken together, 
these could lead to minimum absorption coefficients be- 
low 10~ 8 cm^ 1 per resolution element in 5 s acquisition 
times. Finally, with appropriate change of the DFG and 
upconversion nonlinear crystals along with the CW up- 
conversion pump laser, the system could span much fur- 
ther into the infrared using the same Yb:fiber femtosec- 
ond laser and near infrared VIPA detection. The com- 
bined sensitivity, resolution, wavelength flexibility, band- 
width, and absence of moving parts make this upconver- 
sion femtosecond laser approach a viable and promising 
option for rapid detection of trace gases or separation of 
complex mixtures in a growing field of applications. 
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